clear all;
close all;
x = [0:0.01:10];
x = x+0.05;
z = ((log(x+0.01) - x.^0.95 + 10)./100 -0.087).*300
kir2 = z(find(z>0));
figure; hold on;
plot( kir2);


% y = p1*x^3 + p2*x^2 +
%       p3*x + p4 
% 
% Coefficients:
%   p1 = 6.5808e-07
%   p2 = -0.00028218
%   p3 = 0.029097
%   p4 = 0.067458
% 
% Norm of residuals = 
%      0.27552